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First we present a model independent operator analysis relevant to such decays, then we 
explicitly compute the LFV operator coefficients [and {g^ — 2)] in a general unconstrained 
MSSM framework, allowing slepton mass matrices to have large /i-r entries. We system- 
^ ■ atically study the role and the interplay of dipole and non-dipole operators, showing how 

^ . the rates and the mutual correlations of those LFV decays change in different regions of 

the MSSM parameter space. Values of the LFV branching ratios in the experimentally 
interesting range 10"^ - 10"^ can be achieved. For at least two MSSM Higgs bosons, the 
branching ratio of the LFV decay into fir can reach values of order 10~^. 
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1 Introduction 



The search for lepton flavour violating processes is an important window into physics beyond 
the Standard Model. In the Minimal Super symmetric extension of the Standard Model (MSSM) 
several such processes can have non-negligible rates, provided the mass matrices of the leptons 
and of the sleptons are not aligned. The experimental bounds strongly constrain the amount of 
misalignment in the e-yU sector, whilst the constraints are weaker for the fi-r or e-r sector (see e.g. 
[1,2]). In this work we thoroughly explore the case of large /i-r lepton flavour violation (LFV) in 
a general unconstrained MSSM framework, assuming neghgible LFV in the other sectors. The 
case of large //-r LFV is partly motivated also by the observed atmospheric neutrino anomaly, 
which can be interpreted in terms of — > fr oscillations with maximal mixing angle [3]. The 
prototype and most studied LFV process involving /i and r is the radiative decay r /i7, 
which proceeds via dipole operators. Besides this decay, we also consider other interesting LFV 
processes, such as r — > /lee, r —>■ /x/x/x, r —>■ /ip, r —>■ /itt, t —>■ /ir], t — > pr]', Z —>■ /it, which 
may have different degrees of correlation with t ^ /I'j. The present experimental bounds on 
all these decays are: 
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(Wc arc not aware of experimental bounds on r ^ PV' ■) the MSSM with large fi-r flavour 
violation, a combined analysis of all such processes is very instructive, especially in view of the 
future sensitivity on those branching ratios, which may reach 10~^ — 10~^ [9, 10]. In particular, 
it is interesting to study the interplay between dipole and non-dipole operators. In addition, 
we consider the LFV decays of MSSM Higgs bosons into fir [11], which are correlated to some 
of the above processes and can give testable signatures at future colliders. 

The outline of this paper is as follows. In Section 2 we introduce the relevant effective 
operators and compute the branching ratios of the LFV decays. In Section 3 we focus on the 
MSSM, exhibit the sources of LFV in the left and right slepton mass matrices and discuss the 
computation of the operator coefficients [including the Higgs- /x-r and {g^ — 2) ones], whose 
analytical expressions are displayed in Appendix. In Section 4 we give a numerical discussion 
of the branching ratios in the case of large LFV in the left sector. In Section 5 we present an 
analogous discussion for large LFV in the right sector. Finally, in Section 6 we summarize our 
results. 
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2 Effective operators and branching ratios 



This section is devoted to the 'model-independent' calculation of the LFV branching ratios. 
Namely, first we parametrize the basic effective operators (Section 2.1), then construct appro- 
priate effective lagrangians and compute the branching ratios in terms of the coefficients of the 
effective operators (Section 2.2). We also discuss correlations among different processes in cases 
of single-operator-dominance (Section 2.3). 



2.1 Parametrization of LFV effective operators 

Our first step, in the derivation of the LFV branching ratios, is the parametrization of the 
LFV operators that contain a muon, a tau and either a Z boson, or a photon, or an additional 
/-fermion pair. The leading contributions to these operators arise from o? = 6 SUi^w x 
C/(l)y-invariant operators [12], possibly with Higgs fields set at their vacuum expectation values 
(VEVs). It is useful to keep this in mind, although we will not write exphcitly the operators 
in the unbroken phase (except for a few examples). We also recall that the MSSM contains 
two Higgs doublets Hi^H^-, whose VEVs define the ratio tan/? = {H^jl^H'^. We postpone 
the discussion of o? = 4 Higgs-muon-tau effective operators to Section 3.2. The parametrization 
given below assumes fermions to be on-shell, whereas gauge bosons could also be off-shell. 
We keep the tau mass m,- at the leading order and neglect m^^m^ as well as the light-quark 
masses. We adopt two-component spinor notation, so for example /i and r [jiF and f '^) are the 
left-handed (right-handed) components of the muon and tau fields, respectively^. Sometimes 
{e.g. in Figs. 1, 2, 4) symbols like iJ>,T,f will generically refer to the particles, not to specific 
chirahties. Finally, it is understood that the coefficients of all the LFV operators below should 
carry a flavour subscript /it, which we omit for brevity. 

T/iZ effective operators 



We distinguish the following operators: 

gzml [aI Jta'^T + A| //V^r^ + h.c] Z^, (8) 

gz [Ci Jla'^T + C| /xV^r + h.c] DZ^, (9) 
gzrrir [iOf JIa^"'r' + iD^ ifa^^'r + h.c] Z^^, (10) 

where gz = \/ g^ + g'^ {g,g' being the SU{2)w and U{1)y coupling constants, respectively) 
and Zfj^i, — dfjtZi, — di,Z^. The operators in (8) are chirality conserving (monopole) and have 
no derivatives, so they originate from SU{2)w x C/(l)y-invariant operators with at least two 

^For instance, in terms of the four-component spinors Vj = (m 'P")^ V't = (''' ^'^)) the bilinears Jla'^r and 

ifa^r'^ correspond to tp Pl'^I-'t and i)^^YPR^'T [Pl,r — ^(1 T 7^)] respectively. We take cr'' = (l,cr), (T^ = 
(1,-<t), a^"" = l{a''a''-a''ai'),a'"' = \{a>' a" - a"" a^") and c/^^ = diag(-(-l, -1, -1, -1), where = CT^ ct^) 
are the 2x2 PauU matrices. 
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Higgs fields, which refiects in the factor m| we have extracted out. For instance, the leading 
contributions to (8) come from operators such as: 

{L^a^L,) (hIiD^H,) , (hUSD^H^) , {^I'a^r) {hUd^H,) , (11) 

and analogous ones with hIHi replaced by H2H2 or if2-f^i(+h.c.), where L^, Lr are the SU {2)w 
lepton doublets. The operators in (9) are also chirality conserving (monopole) and have two 
derivatives, so they originate from SU{2)w x ?7(l)y-invariant operators with no Higgs fields 
at the leading order. The operators in (10) are chirality flipping (dipole) and come from 
SU{2)w X C/(l)y-invariant operators with at least one Higgs field, hence we have pulled out 
a factor m^. In the case of virtual Z, only the A^-operators need to be considered since the 
remaining ones (and other ones not listed above) give suppressed contributions. In the case of 
real Z, the A^- and C^-operators give contributions of the same order^ (notice that □ — > —m\). 
The D'^-operators give comparable contributions if the m,- suppression is compensated by a 
large tan /3 factor, induced by the Higgs field H2 (notice that 50 x m-r ~ mz) ■ 

effective operators 



We distinguish the following operators: 

e [CI Jla^'r + n^a^r'' + h.c] OA^, (12) 

errir [iDl Jla^''r^ + iDl n''a^"'r + h.c] F^^. (13) 

The operators in eqs. (12), (13) are analogous to the and operators, respectively. We 
recall that the Z and photon operators have a common origin in terms of SU{2)w x C/(l)y- 
invariant operators. However, the A^-operators (8) have no photon counterpart because of 
the unbroken [/(l)em symmetry. In the case of virtual photon, both the monopole and dipole 
operators contribute, whereas in the case of real photon only the latter ones do. 

Tfiff effective operators 

We parametrize four-fermion operators as follows^: 

^ [iW'r) (B-t l<y,f + Bi- faj") + (/.V^r^) (^^ Jc7,f + B^^ fa.T) + h.c] . (14) 
/ 

Here we have retained chirality conserving operators only and neglected chirality flipping ones, 
e.g. {n'^T){f f^), because the latter are expected to suffer from a double chiral suppression 
and so would give subleading contributions to our LFV processes (in particular this holds in 
the MSSM, even for large tan/9). We emphasize that eq. (14) does not yet include either 
the contributions induced by Z and photon exchange, to be added in the next section, or the 
chirality flipping operators induced by Higgs exchange, which are relevant for large tan (5 and 
will be discussed in Section 3.2. 



^Sometimes this feature has been overlooked in the literature, by assuming the chirality conserving operators 

to be dominated by the zcro-momcntiim component, i.e. in our notation (sec e.g. [13]). 

■^Another discussion on the phenomenological implications of fi-T LFV four-fermion operators appeared in 
[14]. 
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2.2 Effective lagrangians and branching ratios 

All the coefficients of the operators hsted above (A^, , , , D'^ , B^) have mass dimension 
—2, and encode information on the underlying physics. Such coefficients are the building blocks 
for the amplitudes of our LFV processes, as schematically shown in Fig. 1. Now we present for 
each LFV decay the relevant lagrangian and the corresponding branching ratio. 








/i 
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Figure 1 : The different contributions to r — > fi'y, Z ^ and r — > fif f decays. 



The effective lagrangian for the Z ^ jir decay is easily obtained from eqs. (8-10): 



where 

The branching ratio reads as: 



pZ — _ fiZ 



I pZ|2 I I pZ|2 I 1 



m. 



-r T^Z 



mz 



+ 2 



mz 



+ h.c., 
(15) 

(16) 



BR{Z tr) (17) 
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where c = (1/4 - + 2s^y^ ~ 7.9, BR(Z ~ 3.4% and lepton masses have been 

neglected in the kinematics^. 

r ^ //7 



In this case the expression of the effective lagrangian is already in eq. (13). The related 
branching ratio is: 

BR{t- ^ /.-7) = ^ \\Dl\' + \Dl\'] BR{t- /.-P,^^ (18) 



F 



where a = e^/(47r) is the fine-structure constant, Gi;' is the Fermi constant and BR{t 
r — ^ fii'^i", T — ^ fip, T — ^ fiP {P = IT, r], T]') 



The four-fermion effective lagrangian relevant for these decays is obtained by combining the 
contributions in eq. (14) with those induced by Z and photon exchange (see also Fig. 1). Other 
contributions, induced by Higgs boson exchange, will be discussed in Section 3.2. We have: 

O = E[(7^^''^)(^/'7^M/ + ^/Vvj') + (/^w)(Fiv^^/ + i^iVvJ7 
/ 

- 2e^ {Dl TI-a^^T^ + Dl ^l'^a^•^r) ^ ^ Qf(faJ + PaJ") + h.c. (19) 

/ 

where 

FtiR) = Bi\R) + iTl-Qfsl,)glAl^^^ + Qfe'Cl^^^ (20) 

HtR) = Bl'lR)-Qf^w9lAHR) + Qfe'Cl^Ry (21) 

For the sake of brevity, wc will often refer to the operators in the first line of eq. (19) as 
'monopolc' ones, to be contrasted with the 'dipolc' ones in the second line. 

The effective lagrangian (19) can be applied either to leptonic transitions, such as the decays 
T" — > ^'e'^e" and r~ or to semileptonic transitions involving hadrons, after 

taking the appropriate matrix elements of quark bilinears. In the processes involving a neutral 
pseudoscalar meson P {P = tt, 77,77') the relevant matrix elements are (0| J5'"(0)|P(p)) = ifpp^ 
(a = 0,3,8), where fp are the decay constants and Jg*" are the axial-vector currents^. Among 
the neutral vector mesons, we have selected the p meson and so the relevant matrix element is 
(0| J^'^(O)Ip) = Kpm?pe^, where J^^ is the third component of the vector current, Kp ~ 0.2 and 
mp, are the p mass and polarization vector, respectively. Our results can be easily extended 
to tlie otlier A'cctor rnesoiis. 6 and u. 



"^Notice that the factors TOi- in eq. (17) come from the D^-operators, which may contain a (compensating) 
large tan/3 factor, as mentioned above. 

^In terms of four- component spinors V'J — {q '<f ) {q = u, d, s), collected in a triplet = {tp'^ Vj V'f )> the 
vector and axial- vector currents read as J'^" = ^'7'*^^' and J^" = ^'7'*7^^^' (a = 0, 1, ... ,8), respectively. 

The neutral currents contain the matrices A° = y^diag(l, 1, 1), = diag(l, —1, 0), = ^diag(l, 1, —2). 
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Iiee, T — > njiji 



The branching ratios are directly computed from eq. (19): 

_ ^ 1 



BR{t- 



H e e 



8Gl 



eL\2 



\Fl-\' + \F, 



R I 



+4e^Re (i^KF^ + Fl^) + Dl{F'^- + 



+8e^(|D2|^ + |D^|^) log^-3 



BR{t- 



(22) 



BR(r' 



jJL jJL jJi' 



+4e2Re (i?I(2Fr + ) + mK"^ + 2^^")) 
+8e^(|Dlp + |Djr)flo, 



mf 



5i?(T- ^ ^-y^^ur). (23) 



The electron mass rUe and the muon mass regularize an infrared singularity in phase 
space in T ^ ^ee and r — > respectively. 

The effective lagrangian for r ^ up is also obtained from eq. (19): 



'"TUP ~ '^P 



where 



m 



{Fljia^T + F^p'a^'r') + 2e^mr {Dlpa^^^r + L'^/xV^'r) + h.c. 

(24) 



Fl{r) 



1 1 

2^ L{R) ^ L{R)) ^ 2^^(^) ^L{R)) 



i^LiR) ^L{R), 



(Blfn) - BtiR)) + (1 - 2^^)5|Af(^) + 2e2C2(^) 



(25) 



The related branching ratio is: 



+ 



e^{2 + x 



\F[\' + \K\' 



6e^ 



1 + 2t 



-Re {dIF£ + DIF^) 



{\Dl?+\Dl 



,7|2> 



Xp(l + 2xp) 

where = cos^ 6'c c± 0.95, Xp = TJip/m'^ ~ 0.19 and BR{t~ i^tP~) — 25%. 



(26) 
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T — > /XTT 

From eq. (19) we also derive the effective lagrangian for r — > /xtt. It does not receive 
contributions from photon exchange, and reads as: 



^t. = U {FlW'r + F^i^^a'^r) + h.c. 



(27) 



where f^ = fl'^ 92 MeV (we neglect /°, /^) and 



jjiTr 



1 



p"!, \ _( TP'>J'R IT" 



L(R) -^L{R) 



2 aZ 
L{R) 



The branching ratio is: 



(28) 



BR(t- //-tt") = 
where BR{t~ — > Ur'^") ~ 11%. 



^GWc 



(29) 



• T 



Each of the two mesons 77 and rj' has both octet and singlet SU (3)fiav components. The 
relevant decay constants are /p, where P = rj,rj' and a = 8, (we neglect fp) [15, 16]. 
For the ratios fp/f-n we take the following numerical values [16]: ~ l-2/7r, /° ~ 
0.2A, /8 ~ -0.45A, ~ 1.15A. 

The effective lagrangian for r — > //Ty, which does not receive contributions from photon 
exchange, is derived from eq. (19). It reads as: 



{f,Ff + f^Ff) W'r + {f^F^ + /.V^r] + h.c. (30) 



where 



^L{R) 



^L{R) 



1 



K^LiR) 



Fr'Tr,, -2F, 



sl 



1 



( 'P^B _|_ P"-R 9 ^ 

L(iJ) L(R)) 2\/3 ■* ■'^(■R) L{R)/ 



2^/3 
1 



+ ^L{R) - '^^lIr)) - i^ltR) + Bl{R) - ^^L{R)) + 9z^ 



2 aZ 
L{R) 



(31) 



1 

76 



ul 
L{R) 



+ Bi\p) + B 



L(_R); 



+ 5 



• (32) 



The branching ratio can be expressed as: 



BR{t- iJ-r)) 



f8 fO 
•'11 Tpri,8 , Jv 771)7,0 



+ 



+ {L^R) 



[1 - Xr,f BR{t- i/^TT") 



(33) 
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where Xn = m^/m^ ~ 9.5 x 10 ^ and the ratio m^/ml has been neglected. For the decay 
T — > ^T]' , one has to make the replacement 77 77' in eqs. (30), (31), (32) and (33). We 
anticipate that the equalities F^^j^) = -^z!(^)) ^l(r) ~ ^l{r)^ which follow from eq. (19), 
will not hold in the presence of sizeable Higgs-mediated contributions [see Section 3.2]. 



2.3 Correlations 

In the previous subsection we have seen that, on the one hand, each operator can contribute to 
different processes, hence correlations exist. On the other hand, the correlation pattern is not 
trivial because each process can get contributions from different operators, so making general 
predictions is not straightforward. Nevertheless, it is interesting to see what happens when 
several processes are dominated by the same operator. Such a situation may or may not be 
realized in a specific underlying model or in a specific portion of the model parameter space. 

D'^ -dominance. If photon dipole contributions are the dominant ones in the decays r — > 
//ee, T — > /////X and r — > //p, then such processes can directly be compared to r — > 



BR{t —>■ 11 e^e ) a ( m 



2 



— log^-3 -10-^ (34) 



BR(t — >• /i 7) Stt \ m 

^ ^ ^ flog^ - 2.2 X 10- (35) 

BR{t- fi-j) 37r V 4 y ^ ^ 

BR(r- ^ fi-p^) ^ a (2 + x,) BRjr- ^ u^p-) ^ o 3 ^p-a .ggx 
BR{t- /i-7) UttcI Xp{l + 2xp) BR{t- i^-p^Ur) ' 

In particular, the present bound (1) on r — > //7 would imply the following bounds: 

BR{t~ /i-e+e") < 3 X 10"^ (37) 

BR{t- /i-/i+/i-) < 7 X 10"^°, (38) 

BR{t- 11- p^) < 8 X 10"^°. (39) 



-dominance. Suppose now that ^'''-contributions are suppressed and that monopole con- 
tributions to T ^ //ee, T — > pLpLpL, T ^ /ip are dominated by C'''-operators. Then: 

BR{t- /x-p°) ~ BR{t- p^'l^^p^') ^ 1.5 X BR{t- n-e^e-). (40) 

A'^ -dominance. Suppose again that dipole contributions are suppressed and that monopole 
contributions to r ^ jiee, r — » fififi, r — > /ip, r — > /itt, t — > /ir), r — > /irj' and Z ^ /it 
are dominated by A'^-operators. Then: 

BR{Z fx+T~) ~ 3 X BR{t- /i"e+e"), (41) 
BR{t- pTp^) ~ 1.8 X BR{t- /i"e+e"), (42) 
BR{t- //-7r°) ~ 2.7 X BR{t- l^'e+e'), (43) 
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BR{t- ji-ri) ~ 0.8 X BR{t- A*~e+e"), (44) 
BR{t- fx~r]') ~ 0.7 X BR{t- f^'e+e'), (45) 

BR{t- ^ fx-ix^fx-) ~ Q'^) X 5i?(r" ^ /i"e+e"). (46) 

The result for r — > /i/iii depends on the relative amount of Af and contributions: 
the upper (lower) estimate refers to the case of pure Af (A^) dominance, with A| = 

(^! = o). 

B^ -dominance. Finally, consider the case in which all contributions are suppressed except 
for those S-^-induced. Then correlations are in general weakened or lost, since differ- 
ent B'^ coefficients appear in different processes. In some specific case, however, simple 
correlations may emerge. For example, consider the limit of fiavour universality in the 
down-quark sector, i.e. E^^^) = 5^^^^ and = 5^^^) (which holds in the MSSM 

with down-squark mass universality). T?hen we obtain BR{t — > /xtt) ~ 3 x BR{t — > ^rj) 
by neglecting /° in eq. (33). Another example of correlation is that between t ^ /ip and 

r — > /iTT m the absence of LFV in the left sector (E^^'""^ = 5/*'^' =0). In this case, m 
the S'C/(2)vi/-symmetric limit (-B^^ = B'j^), both decays only depend on the combination 
B^" - hence BR{t fip) - 2.3 x BR{t pir). 



3 The MSSM framework 

In the above sections we have given a model-independent description of effective operators that 
contribute to our selected set of LFV processes. Here our purpose is to compute the coefficients 
of such operators in the framework of the MSSM. The source of LFV in the MSSM is the 
potential misalignment between the lepton and slepton mass matrices. More precisely, in the 
superfield basis in which the charged lepton mass matrix is diagonal, the sources of LFV are the 
off-diagonal entries of the soft-breaking matrices M\,M.\ and A. If we restrict our attention 
to the second/third generation, the latter matrices have the form (see also Appendix A.l): 



Ml 



Ml=\ ' I- (4^) 

^iV ^Irrl yrriR^CT rhRrrJ V ^r^^^ KA^ 



Notice that we have parametrized the entries of A by extracting a suitable Yukawa coupling 
(either or hr). We will use the symbols (LFV)l and (LFV)r to characterize LFV in the left 
and right sectors, respectively: 

(LFV)^: mi^.^0 and/or ^ 0, (48) 

(LFV)k: ml^^j^O and/or ^ 0. (49) 

In the case of pure (LFV)l [(LFV)ij] the parameter A^^ (^^r) expected to be suppressed by 
h/x/hr with respect to Aj^^ {^^t)- 
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In Section 3.1 we describe some general features of the MSSM contributions to the operator 
coefficients hsted in Section 2.1. We also present the Higgs-^u-r effective operators and discuss 
their impact on r ^ r fin, r fir], r firf' (Section 3.2). Finally, we discuss the 
effect of LFV on the muon anomalous magnetic moment (Section 3.3). The explicit results of 
the calculations for the operator coefficients are presented in Appendices A.2-A.10. 

3.1 General aspects of operator contributions 

The leading contributions to the coefficients A^, , , , D^, arise from one-loop dia- 
grams that involve the exchange of gauginos, Higgsinos and sleptons (and also / sfermions 
in the case of B^). These diagrams are shown schematically in Fig. 2 and in more detail in 
Appendices A.3-A.8, where the explicit analytical results for the coefficients are also presented. 
The relevant parameters involved are those in eq. (47) (or equivalent ones, see below), the 
SU{2)w X U{1)y gaugino masses (M2, Mi), the fi parameter^, other sfermion masses and tan/5. 
We neglect contributions to the LFV operator coefficients coming from other sources^. 

In the diagrams in Fig. 2, the flavour transition occurs along the slepton line. The generic 
behaviour of the related operator coefficients at the linear order in the LFV parameters m^^, 
(the so-called mass-insertion approximation [18]) goes as follows: 



where the constant g^, can be either g or g' and M5 is some effective sparticle mass. These 
relations show that generically all the coefficients have the same order of magnitude (except 
for a possible tan /9-enhancement in the dipole coefficients), are suppressed by M| and are 
proportional to the relative amount of LFV. Nevertheless, detailed computations arc, of course, 
necessary to single out the specific dependence of each coefficient on the MSSM parameters. 
Upon addressing such a computation one is concerned, in particular, with the treatment of LFV. 
Since LFV in the second/third generation may not be a small effect, we choose to go beyond 
the mass-insertion approximation. Specifically, we diagonalize the mass matrices M.\,M.\ and 
work with their mass eigcnstatcs, i.e. our results are expressed in terms of the corresponding 
eigenvalues fh\^,rh?^^ {a = 2,3) and mixing angles 9l,0r. In particular, large (LFV)l in 
means m|^^ ~ '"^1^,* ~ ''^\tt-i which in turn implies 6^ ~ C'(l) and {rh\^ ~"^L3)/(^i2 '^'^'lz) ~ 
0{1) (with ?7i|2,m|g of the same order or hierarchical). Analogous relations hold for the case 
of large (LFV)/?. Another aspect we have to deal with is the inclusion of clcctroweak breaking 
effects. We choose to treat such effects at lowest order, i.e. the only Higgs insertions we 
consider are those explicitly depicted in the diagrams in Fig. 2 (and also in Appendix), which 

^We take a superpotential term of the form W = iJ,HiH2 = /x(-ff°if2 — ^2 ), so our sign convention for 
the fj, parameter is opposite to the one commonly used (for this and other conventions see also Appendix A.l). 

^Additional contributions arise, for instance, from the interactions of the MSSM fields with the goldstino 
supermultiplet [17]. Such effects are suppressed, unless the supersymmetry breaking scale a/F is close to the 
electroweak scale. 




(50) 



(51) 
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X ~ X 

X\H\\ 



+ 



~ X 
H \ A 



X 



+ 



yl^(c). 

A 



X 




A 



+ 



~ X 
H \ A 



+ 



A 



X 




A 




A A 



Figure 2: Topologies of the MSSM diagrams contributing to the coefficients , , {V = 7, Z) and . 
Internal solid lines denote gauginos (A) or higgsinos {H). Dashed lines denote sleptons (or other sfermions in 
the box diagrams). Dotted lines with a cross denote Higgs insertions. 



also corresponds to take only the d = 6 operators in the SU{2)w x t/(l)y-unbroken phase. 
We believe that working in this approximation is more transparent and appealing from the 
theoretical point of view, although this may imply some loss of numerical accuracy when M| 
is very close to m|. 

A few additional comments are in order about some operator coefficients. In particular, we 
recall that is associated to operators like those in eq. (11), where the field is contained in 
a covariant derivative acting on a Higgs field {D^H). Thus the operator coefficient can also be 
extracted from the part of the operator containing an ordinary derivative ((9^ if), that is from 
diagrams with momentum dependent Higgs lines and no Z lines. This method, which we have 
used, has the advantage that one need not compute LFV diagrams corresponding to wave- 
function renormalization of lepton fields, which are associated to operators like {L^^ DyLT-) 
with or without Higgs fields attached. Also notice that receives contributions of three 
different types: 

aZ _ AZ{a) .Z{b) .Z{c) /c-oN 

The terms Af^^^ , A^^ are parametrically suppressed by / with respect to A^^^ , hence those 
contributions can only be relevant for large tan /3. Also the dipole coefficients (V — 7, Z) 
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receive contributions of three different types: 



Dln = DV;S + Dl^^ + DV:l (53) 

Notice that the diagrams contributing to D^^^ can only contain the Higgs field Hi, which 
comes from the r equation of motion, whereas either Hi or H2 can appear in the diagrams 
contributing to D]^^^ , D]^''^ . As a consequence only the latter coefficients can receive a tan/3- 
enhancement. In the case of it is sufficient to compute these tan /3-enhanced terms, as they 
are the only ones which can give contributions comparable to the monopole ones in the Z — > //r 
decay [see eq. (17)]. Notice that tan /3-enhanced terms have a more dramatic effect in the case 
of D'^ [19], since they generically make the dipole operators dominate over the monopole ones 
in decays such as r ^ /lee, r — > /i/iii and r /ip. 



3.2 Higgs- muon-tau effective interactions 

At large tan/9, another class of LFV interactions is relevant, namely those between a //-r pair 
and Higgs bosons. For our purposes it is sufficient to focus on the leading effects, which arise 
from these dimension-four operators: 

-hriS,lHl*T^li - hrARH^yr + h.c. (54) 

The corresponding diagrams are sketched in Fig. 3 and shown in more detail in Appendix A. 9, 



A(&): 



-"2 
H \ A 



AW: 



A 



HI 



Figure 3: Topologies of the MSSM diagrams contributing to the A coefficients. Dashed hnes denote sleptons. 

where the explicit expressions of the dimensionless coefficients Al^r = A^^^^ + A^^^ are given^. 
In the mass-eigenstate basis for both leptons and Higgs bosons, the LFV couplings read as: 

where A is the physical CP-odd Higgs field, a is the mixing angle in the CP-even Higgs sector 
[V2Rc{H^ - (H^)) = CaH-Sch, \/2Rc(iJ^ - (iJ^)) = SaH + Cah], and a short-hand notation is 
used for c^ = cos^, s^ = sin^ (^ — a, P, P — a). The effective couplings (55) directly contribute 

*We should add that most of the properties discussed in this section do not rely on the specific MSSM origin 
of Al,r. The discussion can also apply, more generally, to two-Higgs-doublet models in which lepton masses 
are mainly generated through Yukawa couplings to Hi, with smaller contributions from H2. 
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to LFV decays of the neutral Higgs bosons, A,H,h^ fir [11], and also induce Higgs mediated 
contributions to several LFV r decays, such as r — * 3/i [20, 11], r ^ firj [21], r fit]' and 
r — s> /ivr. The A-contributions to all such processes are displayed in Fig. 4. Here we discuss 
each of them. 




Figure 4: A-contributions to the Higgs boson decays A,H,h /ur and to the decays t 3/z, t /Utt, 
T — > fXT], T — > ;ur/'. In the last diagram, curly lines denote gluons. 



Higgs^ IIT 

Concerning the Higgs boson decays, we have 

BR{A i^+T-) = tan^ (3 {\Al\^ + \Ar\^)BR{A r+r") , (56) 

where we have approximated 1/c^ ~ tan^ (3 since non- negligible effects can only arise in 
the large tan/? limit. If A is replaced with H [or h] in eq. (56), the r.h.s. should also 
be multiplied by a factor (c^_Q/sa)^ [or (s^_a/ca)^]- These LFV decays and the related 
phenomenology have been extensively investigated in [11] (for other studies in two-Higgs 
doublet models see e.g. [22]). We recall that BR{A — > fxr) can reach values of order 10~^. 
The same holds for the 'non-standard' CP-even Higgs boson (either H or h, depending 
on itia)- 



Consider now the implications of virtual Higgs exchange for r — ^ yUyUyU. The effective 
lagrangian C^^^ [see eq. (19)] receives the following additional contribution from Higgs 
exchange: 

= -^'^F ^^^^r (C+JI r + (Ar//V + A2//T^) + h.c, (57) 

where the C± coefficients are: 

4 V T^H rn\) 
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The operators proportional to C+ can be Fierz-rearranged and cast in a form already 
exhibited in eq. (19). In other words, these contributions can equivalently be regarded as 
a shift in and Fg^, that is 

5F^^ = ^ m^m,C+A«, 5F^' = ^ m^m,C+A^, (59) 

and the branching ratio (23) is accordingly affected. The operators proportional to 
C- have a different chiral structure and give an extra contribution (m^m^C_ )^(|AlP + 
|Ai?p)/32 to BR{t 3fi)/BR{T fiuu). Notice, however, that the Higgs-mediated 
contributions are potentially relevant only for large tan/?, so that the lepton mass sup- 
pression is overcome. More precisely, the truly enhanced operators are those proportional 
to C+, which becomes C+ ~ 2tan^ f3/m\. In the hmit in which the Higgs-mediated con- 
tribution dominates over the other ones (in particular, the dipole one), the decay r — > 3/x 
is directly correlated to the decay A ^ fir [11]: 



: , , ,/tan/5V nOO GeV\^ BR(A^ fi+r-) , , 



T^^P {P = n,n,n') 

Now consider the implications of virtual Higgs exchange for the decays r — > /iP, where 
P is a neutral pseudoscalar meson (P = tt, 77,77'). Since we assume CP conservation in 
the Higgs sector, only the exchange of the A Higgs boson is relevant. Moreover, in the 
large tan (3 limit, only the A couplings to down- type quarks are important. These can be 
written as: 

-i(\/2G'F)^/^tan/3 Ai^amad^d + CsrUsS^s + Cbmbb^b) + h.c. (61) 

The parameters ^d,S,s,^b are equal to one at tree level, but can significantly deviate from 
this value because of higher order corrections proportional to tan/3 [23, 24], generated by 
integrating out superpartners^. At energies below the bottom mass, the 6-quark can be 
integrated out, so in eq. (61) the bilinear —imi,b'^b+h.c. is effectively replaced by the gluon 

2 

operator Q = gf^e'''''"^G^^G'^^, where Qs and G^^ are the SU{3)c coupling constant and 
field strength, respectively (see Fig. 4). The effective lagrangian due to A boson exchange, 
relevant for r —>■ iiP (P = tt, 77, 77'), reads: 

SC'^^ = iV2Gf 7n,^^ kdmdfa + C.m.jf + i.Vt] {Arij,'t + Ali^r') + h.c, (62) 

777.^ 



^In the limit of quark flavour conservation, each {q = d, s, b) has the form = (1 + Ag ta.n/3)~^, where 
appears in the loop-generated term —hqAqH2*q'^q + h.c. Wc recall that the leading contribution to Ag arises 
from a gluino-squark loop and reads Ag ~ — ^/iM3/3(M|, m?^, m?^). A Higgsino-stop loop also contributes 
to A{, through a term —^iiAi;I^{ii^,m^ ,m3 ). 
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where we have defined the quark pseudoscalar densities^° = i{q — q'^q), with q = d,s, 
and approximated 1/c^ — tan^/3. The matrix elements {0\jg\P) and (0|r2|P) can be de- 
termined along the lines of [25], i.e. by taking the matrix elements of d^J^^ , dij,J^^, d^J^^ 
(the latter divergence is anomalous and contains the term —\/QVl). We find: 

(0|m,j,^|7r) = (OKjfk) = (0|l^k) = -^([^)Am^ (63) 

(0Kjt|,7) = -^/X, (0|Q|,7> = (^/^ + /;) mj, (64) 



(65) 



where 2; = mu/md- In deriving the matrix elements with 77,77' we have neglected mu,md- 
All these results translate into additional contributions to the efi^ective meson lagrangians 
Cfjj^^, Cf^j^^, Cf^j^^, in Section 2.2. Such contributions can be interpreted as shifts in F^j^-^, 

7-ir),8 T-i»7,0 pV.S rpv',0 4-V,„^- ;„ 
-^L(il)' -^L(fl)' ^^^'^^ 

5F,,^, = -V2G.m,(^e.Y^ + ^3^J— ^A,(^), (66) 

^^i.) = mje.^ A,(«) , (68) 

and SF2(^^y '^-^l(^) obtained from ^-^^'(r)' ^^l{r) substitution rj ^ r]'. We 

have performed a cross-check by using the approach of chiral perturbation theory and 
obtained results that are consistent with those above. 

In the limit in which the processes r — > 3// and r — > //jy are both dominated by Higgs- 
exchange, these decays are related as: 



BR{t n fi+jj, ) \m^m 




(69) 



For ~ ~ 1 this ratio is about 5, but it could also be a few times larger or smaller 
than that, depending on the actual values of ^si^b- Our result confirms the relevance of 
Higgs-exchange in r — >^ /xry, previously emphasized in [21]. In the latter paper, however, 
neither the contribution of the (bottom-loop induced) gluon operator Q nor the factors 
were included. Moreover, if these effects were disregarded, the ratio (69) would be 3 
times smaller than the one found in [21] (notice that F^ in [21] corresponds to -\/2/®). 



^°In four-component notation the pseudoscalar densities read as = itjjg'y^ipg with ip'^ = [q q"). 



15 



Finally, let us compare r 
domination: 



/jiT)' and r — > //TT with r — > //Ty in the limit of Higgs-exchange 



BR{t- ii-T)') 



BR{r- 
BR{t- 



H rj) 



H TT) 




BR{t- ji-q) 




II 3 V /Cs , 1 
?6 ^ 3 + 3/8 



1 I in _|_ 6U-2^ 

^ I ]_ I \/2 /?7 

^6 ^ 3 ^ 3/8 



(70) 



(71) 



Both ratios are suppressed, although for different reasons. The ratio BR{t — > /x?]') / BR{t — 
/IT)), which seems to be 0{1), is much smaller because the singlet and octet contributions 
tor yU?7' tend to cancel against each other (we recall that f^i/flj/ ~ —0.4). ForC^s/Cb ~ 1) 
the ratio (70) is 6 x 10~^. The ratio BR{t ^tx)/BR{t fir]) is small because it is 
parametrically suppressed by m^/m^ ~ 10~^. The actual numerical value is sensitive to 
the parameter z — rriu/md, and to the ratios ^d/Cb, Cs/Cb- In the case Cd/Cb ~ Cs/Cb ^ 1; 
if we let z vary between 0.2 and 0.7 [7] the ratio (71) ranges from 4 x 10~^ to 10~^. These 
results, combined with the present bound (6) on r — firj, imply that the Higgs mediated 
contribution to BR{t /ir}') and BR{t — > //tt) can reach (9(10~^). 



3.3 Muon anomalous magnetic moment 

The flavour changing dipole operators in eq. (13) have an obvious flavour conserving counter- 
part, which for the muon is 

em^iDl fi^'aP^fi F^^ + h.c. (72) 

The coefficient DJ^ is directly related to the anomalous magnetic moment of the muon: = 
{g^ — 2)/2 — 2mf^ReD^. In the MSSM the coefficient D'J^ receives three types of one-loop 
contributions from superparticle exchange, analogous to those shown in Fig. 2. In turn the (a)- 
type and (b)-type contributions can be further distinguished according to the chirality (either 
L or R) of the sleptons circulating in the loop. In the (c)-type contributions both L and R 
slcptons are simultaneously present. Therefore we can express the superparticle contribution 
to a^, denoted as a^^^^, as the sum of five terms: 

a^- = + a^;l + af, + a% + a%. (73) 

Non- vanishing contributions to a^Jf^^*^ arise even in the absence of LFV, of course [26, 27]. In 
this case the only sleptons involved are = (i/^,/!^) and (Ir, and the contributions to are 
generically proportional to (g'^/167r^) (m^/M|), possibly with an extra tan/3 factor. In the 
case of mixing between second and third slepton generations the no-mixing results need to be 
generalized (see Appendix A. 10). One of the effects of LFV is that (c)-typc contributions have 
extra terms proportional to m^/m^, as already observed in [27, 28]. These contributions are 
potentially the largest ones if (LFV)l and (LFV)r are both large. However, the latter situation 
does not seem very natural if the smallness of mij_/mr is to be explained by an underlying 
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flavour symmetry. In fact, we rather expect LFV to be large in at most one sector (L or R). In 
such a case the apparent m^/m^-enhancement is compensated by suppressing factors, so that 
the size of such extra contributions to a!!^^^^'^ does not exceed that of the other ones. 

We also recall that, at present, it is not clear whether a discrepancy exists between the 
experimental determination of and the Standard Model prediction (see e.g. [29]). Hence 
some caution is needed in deriving constraints on the MSSM parameter space [30] . 
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4 Large (LFV)^: numerical analysis 



In this section we perform a detailed numerical analysis in the case of large (LFV)^, assuming 
vanishing (LFV)r, i.e. fh\^^ = A^^ — 0. We remark that all operator coefficients depend 
on mL„, 9l and gaugino masses. Some coefficients also depend on additional parameters. In 
particular, and D'l^''^ depend on /j, and f3; D]}'^^ depends on /j,, f3, mf^, A^-, A^^; S{ depends 
on rhj. 

The lightest eigenvalue of is conventionally chosen to be m|g (although our formulae in 
Appendix do not depend on such a choice). To enhance (LFV)^, in all our numerical examples 
we will take maximal mixing, 9^ =71/4 {i.e. fh\^^ = m\^^ in and widely split eigenvalues. 

The mass parameters will be varied in such a way that charged sparticle masses be > 100 GeV, 
to respect the LEP constraints [7, 31]. Wc will also check that neutralino production at LEP 
be either kincmatically forbidden or adequately suppressed [32]. Incidentally, we should add 
that we have performed such checks using the full mass eigenvalues (not just those at zero-th 
order in electroweak breaking). For simplicity, we take a common (soft- breaking) mass fhe for 
selectrons of both chiralities, as well as a common mass ihq for first and second generation 
squarks of both chiralities. 



4.1 (LFV)i with small tan/? 

The decay r — > /i7 poses significant constraints on the MSSM parameter space, especially for 
large (LFV)l. Moreover, the decays r — > /xee, r fifJ'/J', r ^ up are often dominated by the 
dipole contribution. One of our purposes is to study to which extent the latter property holds 
in the MSSM parameter space or, more generally, to study the interplay between dipole and 
non-dipole contributions. We start our numerical analysis by considering in more detail the 
contributions to the dipole operator, taking into account that the present bound (1) on r — > /X7 
translates into the bound \Dl\ < 5 x 10~^ GeV~^. 



Of"'^ [10-9 ■ GeV-2] 



[GeV] 
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150 200 
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Figure 5: Dl^^^ isocontours with tan/J = 3, niL^ = 1 TeV, = 7r/4 and /x = 140 GeV (1 TeV) in the left 
(right) panel. The solid and dashed lines refer to Mi = 100 GeV and Mi = —100 GeV, respectively. 



The coefficient l-D^ | is at most 10 ^ GeV ^, so it fulffis the bound. In contrast, D]} 
and D'l^'^'' may separately exceed the bound in some regions of the parameter space, especially 
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for large tan i3. In this section we focus on small tan (3. Incidentally, notice that monopole 
coefficients do not have a strong tan/? dependence. In Fig. 5 we show the contours of the 
coefficient D"l^'^^ in the plane {M2,fnL^) with tan/3 = ?>, 6^ = ti / A, rhi,^ = 1 TeV, |Mi| = 
100 GeV and = 140 GeV (1 TeV) in the left (right) panel. For the range of M2 and 771^,3 
shown in the figure, D]^'^^ is well above the bound for small /j, whilst it is of the same order 
of the bound for large ji. Also D^^^^ ^"^^ exceed the bound, depending on the 

range of the extra parameters mf^. At, A^^. Even if D']^^^ and D']^'^^ should separately exceed 
the bound, however, mutual cancellations could bring the total dipole contribution ^2 below 
the bound. To illustrate this, in Fig. 6 we depict the L>2 contours in the plane (74^,mf^) for 



400 - 



300 5; 




Ar [TeV] 



Figure 6: contours for tan/3 = 3, friL^ = 1 TeV, Ol = 7r/4, A^^ = and five choices of 

(mL3,At,M2,Afi) in GeV: (100,140,120,100) (dotted), (100,140,120,-100) (siiort-dasiied) , (100,250,150,100) 
(solid), (100,250,150,-100) (dot-dashed), (200,700,300,-100) (dashed). For each example the two external 
lines correspond to |Z>2I = 5 x 10~^ GeV~^ and the middle one to Dl = 0. 



tan/3 = 3,mL2 — ^ TeV, = and several choices of (m^.^, /i. Mi, M2). For each case we 
show three lines: the two external ones delimit the allowed region (I-D2I ^ 5 x 10^*^ GeV^^), 
the middle one corresponds to full cancellation (Dl = 0). The coefficient Djj''^^ is above the 
bound in the four cases with small while it is below in the example with large /x, hence in 
the latter case the allowed region is wider. In the examples with Mi > 0, large and negative 
At values-*^-^ are required to make the sign of D]^'^'' opposite to that of D]}'^\ We also note that, 
in general, the role of Ar would be shared by Ar and A^^ [see eq. (105)]. 

A remark is in order about the muon anomalous magnetic moment. There are obvious analo- 
gies between the contributions to DJ^ (i.e. to ajf^^''^) and those to D'^, such as a similar diagram- 
matic origin and the dependence on common parameters (see Appendices A. 6 and A. 10). In the 

^^A situation with \ Ar \ » niLs, TOtr could destabilize the scalar potential and induce VEVs for slepton fields. 

To avoid this one can take, for instance, a sufficiently large niA- The inequality m.r\Ar + jitanf]] < \/2'mL3'fhf^ 
is another requirement to avoid tachyonic sleptons (here Ol = 7''/4). This has been verified throughout all our 
numerical examples. Similar considerations apply to the (LFV)ij case in Section 5, with L R. 
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fhfj^ [GeV] 

Figure 7: Branching ratios of LFV decays for tan/J = 3, jfiLs = 1 TcV, niL^ = 100 GcV, 9l = 7r/4, 
A^^ ~ 0. In the upper panel, the remaining parameters are: ^ = 120 GeV, M2 = 150 GeV, Mi = 100 GeV, 
mf„ = = 100 GeV, = = 1 TeV. In the lower panel: ji = 800 GeV, M2 = 120 GeV, Mi = -100 GeV, 
At = 0, rile = 120 GeV, rhq = 1 TeV and m^^ = m?^. The solid horizontal line indicates the present bound 
on BR{t /X7). In the upper example: BR{t firj) 6 x 10~^°, BR{t fxr]') 5 x 10~^°. In the lower 
example: BR{t ^ /^tt) ~ 2 x lO'^^^ ^^(^ ^ ^ ^^^^ ^ ^^z) ^ iQ-12 
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case of large (LFV)l, for instance, we have typically that ajf/^ ~ - 0.02 x (D^^^'^'VCeV-^). 

So if D'l^'^^ saturates the present bound on D'^, then |a|2,| < 10"-'^'^, while if D]}'^^ is ten times 

larger than the bound (which should be fulfilled through cancellations), then |a|f2| ~ 10~^. 
This is of the order of the experimental error on and could either be acceptable on its own 
or play a role in explaining a possible discrepancy between the experimental determination of 

and the Standard Model prediction. The connection between other contributions to ajf^^^ 
and Dl is less direct, since some parameters differently affect each of them. For instance, in 
the {LFV)l case that we are considering here, aj^^^'^ depends also on and A^, which do 
not enter the D^'^oefiicients (at least at leading order). In particular, cancellations among 
different contributions to a^^^^ (analogous to those in Dj) may or may not occur, depending 
on the choice of such parameters. All such features loosen the correlation between aj^^^^ and 
D'^. Since we work in an unconstrained MSSM framework, we content ourselves with verifying 
that the numerical examples of this Section are consistent with values of |a,^^^^^| of order 10~^ 
or smaller. Similar considerations apply to the (LFV)i^ case discussed in Section 5 (where we 
shall not come back to this point). 

After this digression on a^, let us continue the discussion on the LFV operators. It is 
interesting to see how the branching ratios of the different LFV processes behave when the 
Dl dipole contribution varies and the monopole ones are basically fixed. In Fig. 7 we plot 
two such examples, where fi is either small (upper panel) or large (lower panel). The specific 
values of the parameters involved are given in the caption. The parameter that varies on the 
horizontal axis (either At or r/if^) has the only role to induce a variation of the Dl dipole 
contribution. The BRs of Z — > /it, t — > //tt, t /irj, r nrj' are constant^^. In the other 
decays we can observe an interplay between the dipole and monopole contributions, which are 
comparable in magnitude in these examples. For instance, in both examples the pure monopole 
contribution to BR{t fiee) amounts to about 10^^, as we can read in correspondence of the 
dipole cancellation (r — > /X7 'dip'). At the points where t —>■ iij saturates the bound, the pure 
dipole contribution is 3 x 10^*^ [sec eq. (37)] and the combined contribution can reach 5 x 10~^. 
In the case of r — > //p we can also notice a strong interference effect between the dipole and 
the monopole contributions. Consider, for instance, the first example. In correspondence of 
the r /i7-saturation points, BR{t fip) is about 6 x 10~^ (left point) and 1.5 x 10^^ (right 
point). These numbers should be compared with the pure dipole contribution at those points, 
0.8 X 10~^ [eq. (39)], and the pure monopole one, 3 x 10~^. 

The previous examples have shown that monopole operators can play an important role. 
We recall that the monopole coefficients (F^, F^, PP, F'^ , F*?, F'?') receive different contribu- 
tions (from A^, C^, C'', 5-^), which have specific parameter dependences. Notice that all these 
contributions depend on rfiLg and M2, so latter masses should not be too large, if we want to 
avoid a strong suppression of monopole coefficients. The relevance of each contribution, as well 
as the possibility of mutual cancellations, is also influenced by other parameters, such as (in 
A^), the selectron masses (in W^) or the squark masses (in S", B"^, B^). In the flrst example of 



^^In fact, BR{Z — > fir) depends on mf^ through the coefficients A'^^'''''\ D'^^^''^\ Similarly, BR{t — > 
/ztt), BR{t ^1]) and BR{t ^irj') depend on rhf^ through the coefBcients However, for 

small tan j3 these contributions are suppressed. 
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Fig. 7, for instance, the large values of mg and rhq imply suppressed box contributions, while 
small fi implies unsuppressed ^'^-contribution (which combines with those of or C^). In the 
second example, is suppressed by the large fi value, whereas B'^ is unsuppressed because mg 
is small. In the following we shall analyse in more detail the parameter dependence of monopole 
contributions, considering one process at a time. 




100 200 300 400 500 600 700 800 900 1000 
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BR(Z fir) [10-9] BR(Z fir) [IQ-^] 
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Fi gure 8: BR[Z — > /it) for tan/3 = 3, jtilj ~ 1 TcV, 9^ = 7r/4. In the upper panel: m/,., = 100 GcV.mf^ = 
100 GeV, M2 as shown, Mi = 100 (-100) GeV for the sohd (dashed) hnes. In the left (right) lower panel: 
Ij, = 120 (1000) GeV, Ml = 100 (-100) GeV and m^^ = 100 (500) GeV. 



In Fig. 8 BR{Z fir) is studied^^ as a function of fi (upper panel) and in the (M2,m/^._,) 
plane for small fi ( left lower panel) or large /x (right lower panel). The behaviour in the upper 
panel shows the destructive interference of the and contributions. For small // we have 
\Af\ > \Cl\. For large ji, A\ is suppressed and BR{Z — > /it) approaches an asymptotic value 
essentially determined by Cf. The lower panels also show that BR{Z fir) depends on 
{M2,mL^) in a different way for small and large fi. In either case, this BR can reach ~ 10~^. 

In Fig. 9 the monopole contribution to BR{t IJ^l^'lJ^) is drawn as a function of (upper 
panel) and in the {MiifhL^) plane for small 11 (left lower panel) or large // (right lower panel). 
The monopole coefficients _p^^^^' receive contributions of different signs from , C2 and . 
The upper panel shows the changes in the interference pattern, since is substantial for small 
jji and suppressed for large yU, while C2 and B^^'""'^ are constant. The lower panels show how the 
monopole-induced BR{t — > /x/iyu) depends on (M2, 771^3) for small and large yU. We can notice 

^•^For earlier studies on Z ^ jir in the MSSM, see e.g. [33]. A more recent computation, with related 
discussion, appeared in [34] . 
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Figure 9: Monopole contribution to BR{t — > A*/"/") for tan/3 = 3, = 1 TeV, 9l = 7r/4. In the upper 
panel: jtils = 100 GeV,mf^ = mjif, = 100 GeV, M2 as shown, Mi = 100 (-100) GeV for the soHd (dashed) 
Unes. In the left (right) lower panel: ^ = 120 (1000) GeV, Mi = 100 (-100) GeV and m^^ = 100 (500) GeV. 
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yuee) for tan/3 ~ 3, = 1 TeV, 9l = tt/A. In the upper 
= 120 GcV, Ml = 100 (-100) GeV for the solid (dashed) 



Figure 10: Monopole contribution to BR{t - 
panel: niL^ = 100 GeV,TOf„ = 100 GeV, M2 
lines. In the left (right) lower panel: /i = 120 (1000) GeV, Mi = 100 (-100) GeV, wg = 1000 (120) GeV and 
rhff, = 100 (500) GeV. 
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that the maximal monopole contribution to BR{t — > /i/i/x) in these examples is comparable to 
the maximal allowed dipole contribution [eq. (38)], so the combination can be around 10~^ (see 
also Fig. 7). 

In Fig. 10 the monopole contribution to BR{r — > iiee) is plotted in the (/i, mg) plane (upper 
panel) and in the (M2, ifiL^) plane for small and large mg (left lower panel) or large and small 
fhe (right lower panel). The contours in the upper panel reflect the interference effects among the 
different contributions to Fi^'"^\ i.e. Af,Ci and B]^'""^ . The monopole-induced BR{t fiee) 
can be (9(10^^) cither for small fi and large mg or for large /i and small mg. In the former 
case is substantial and box contributions are suppressed, in the latter case the opposite 
situation occurs. In both cases there is the constant contribution. In other regions the BR 
is smaller than 10~^ because of mutual cancellations. The maximal monopole contribution in 
these examples is comparable to the maximal allowed dipole contribution [eq. (37)], and the 
combination can be a few times 10~^ (see also Fig. 7). 



BR(r ^ iJ.p)u-'=o [10-3] BR{r ^ l^p)Dy=o [IQ-^] 




H [GeV] M2 [GeV] 



Figure 11: Monopole contribution to BR{t — > fip) (upper panels) and BR{t — > fin) (lower panels), for 
tan/3 = 3, niL^ = 1 TeV, Ol = 7r/4. In the left panels: ifiL^ = 100 GeV,mfn = 100 GeV, M2 = 120 GeV, 
Ml = 100 (-100) GeV for the solid (dashed) lines. In the right panels: n = 120 GeV, Mi = 100 GeV, 
mg = l TeV and rhf^ = 100 GeV. 



In Fig. 11 wc depict the monopole contribution to BR{t ^p) (upper panels) and BR{t 
/xtt) (lower panels), either in the (iijirhq) plane (left panels) or in the (M2, 772^3) plane (right 
panels). Also for these processes cancellation effects are visible (see left panels). Indeed, box 
contributions interfere destructively with Af and CI- The cancellation regions for t ^ /ip and 
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T jjLTi arc somewhat different, also because C2 only contributes to r — >■ jip. Outside the 
cancellation regions the monopole-induced BR{t — > yup) and BR{t fin) can exceed 10"^. 
We recall that r fip, at variance with r — > pn, also gets a dipole contribution. Notice 
that eq. (26) can be written as BR{t fip) ~ BR^j~,^q^ + BR^p^^^ ± l.h^^ BR^j^^^^BR^p^^^ , 

where BR^j^^q^ {BR^p^Q^) denotes the pure monopole (dipole) contribution. For instance, if we 
combine a value BR^p^^^^ ~ 3 x 10~^ with the maximal allowed BR^p^^^ ~ 0.8 x 10~^ [eq. (39)], 
we obtain BR{t — > jip) ~ 6 x 10~^ or 1.5 x 10~^, depending on the interference sign (see also 
Fig. 7). We have not shown the corresponding examples for BR{t prj) and BR{t — >■ pr]'), 
since they exhibit a similar pattern to that of BR{t pn), though the maximal achievable 
values are somewhat smaller. 



4.2 (LFV)/, with large tan/? and large masses 

Now we discuss the scenario with large tan j3. In this case we have a strong enhancement of the 
jjjib,c} ^jpQ^g coefficients, and Higgs-p-r effective operators become relevant. Two strategies 
can be envisaged to deal with the former issue. In principle, the parameters can be chosen in 
such a way that cancellations occur (as illustrated in the previous section) and keep the total 
Dp dipole coefficient below the bound. However, a significant fine tuning is required, especially 
if M2 and m^g are light. Alternatively, the individual dipole contributions can be kept below 
the bound by taking large mass parameters. This is the option we pursue here. In this case, the 
monopole contributions are suppressed, while the Higgs ones are not (for small m^), because 
the A coefficients are insensitive to the overall mass scale. Hence, at large tan/3 and large 
masses (except for m^i) the following picture emerges. 

i) Z ^ pT is strongly suppressed (even including contributions). 

ii) The decays r — > pee and r ^ pp are generically dipole dominated, so they are correlated 

to T ^ through (34), (36) [and the bounds (37), (39) hold]. 

in) The dominant contributions to BR{t ppp) [eq. (23)] are induced by the dipole term D^p 
and the Higgs-mediated term SFj^^, proportional to A^;, [eq. (59)]. Therefore, r — > ppp 
is correlated to both r ^ pj and A ^ pr [11]. 

iv) The decay t ^ prj is dominated by the Higgs-mediated terms SFp^, SFp^ [eqs. (67), (68)]. 
If r — > ppp is mostly induced by Higgs-exchange, then these processes are correlated as in 
eq. (69). Higgs-exchange is also relevant for the decays r — > //tt and r — > prj' through the 
contribution 6FI [eq. (66)] and 5F/'^, 6F2''° [analogous to eqs. (67), (68)], respectively. 

Let us focus on the processes r — > p^, r — > 3p, A — > pr. The behaviour of these decays 
is illustrated in Fig. 12 in the plane {\p\, IM2I) for tan/3 = 40. Two examples are shown, 
where sign(MiM2) is either positive (upper panel) or negative (lower panel). The relative sign 
between p and M2 is immaterial for large tan f3, while that between Mi and M2 matters because 
it determines the interference between Dj^^^^ and Dl^'^\ as well as that between A^^^ and A^^. 
In the upper panel [sign(MiM2) > 0] the interference is constructive for the Dp components 
and destructive for the A^ ones. Indeed, we can see a cancellation region for r — > 3p and 
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Figure 12: Branching ratios of r ^ ^-y (dashed lines), r ^.fifi (dotted lines) and A ^ jir (solid lines) for 
tan/3 = 40, Ol = 7r/4, ifiL^ = 3 TeV, m^g = 1 TeV, \Mi\ = 5IM2I. In the upper (lower) panel sign(MiM2) 
is positive (negative) and rrif^ = rn^^ = 2.5 TeV (2 TeV). We have also fixed = 100 GeV (which affects 
T 3fi) and BR{A tt) =0.1 (which affects A jir). 



26 



A — > //T [11]. In the lower panel [sign(MiM2) < 0] the opposite situation occurs and we can 
notice a cancellation region for r — * fi'j. In both examples, the contours of r 3/i and A ^ fir 
follow a very similar pattern and are correlated according to eq. (60). This occurs in most of 
the parameter space shown in Fig. 12 where r — > is indeed essentially dominated by the 
Al contribution. Deviations occur in the regions where the D]^ contribution to r ^ 3// is 
comparable to the A^, one (or larger). Regarding the numerical values, in the first example 
BR{A fir) and BR{t 3fi) can reach 5 x 10~^ and 2 x 10~^, respectively. In the second 
example BR{A fir) and BR{t 3/i) can be larger than 10"'^ and 5 x 10^^, respectively. 
Notice that in Fig. 12 we have fixed = 100 GeV. Since the Higgs-mediated contribution 
to BR{t S/jl) scales as l/m\, the values in both examples would be accordingly reduced 
for larger m^. Finally, we recall that BR{t — > /ir)), which can be inferred by eq. (69), is 
generically larger^^ than BR{t 3fi) and can approach the experimental bound (6). In this 
case BR{t — >• /xtt), BR{t jir]') can reach 10~^ [eqs. (71), (70)]. 



5 Large (LFV)i^: numerical analysis 

In this section we perform a numerical analysis in the case of large (LFV)ij, assuming vanishing 
(LFV)/,, i.e. rfi^^^ = A^^ — 0. All operator coefficients depend on niR^, Or and the bino mass 

Ml. Some coefficients also depend on additional parameters. In particular, A^"'^ and D'^^^^ 
depend on ji and /3; D'^^'^^ depends on /5, mf^, At-,A^^; B{^ depends on mj. 

The lightest eigenvalue of Ai'ji is conventionally chosen to be m|.,^ (although our formulae 
in Appendix do not depend on such a choice). To enhance (LFV)^^, in all our numerical 
examples we will take maximal mixing, 9r = ir/i {i.e. = in A4r), and widely split 

eigenvalues. Again the mass parameters will be varied in such a way that charged sparticle 
masses be > 100 GeV. The parameter Mi will be allowed to reach smaller values without 
being in conflict with LEP bounds on neutralino production (e"'"e~ — >■ XiX2(3)) non- 
standard contributions to the invisible Z width {Z XiXi)- We also recall that in the (LFV)/? 
case the parameter M2 does not play a direct role and that we do not impose a speciflc relation 
between Mi and M2, hence small Mi and large M2 can coexist. It is worthwhile to explore 
low Ml values because the monopole coefficients are enhanced. For simplicity, as assumed in 
Section 4, we take a common mass mg for selectrons and a common mass rfiq for first and second 
generation squarks. 



5.1 (LFV)i^ with small tan^^ 

In analogy with the (LFV) l case, first we focus on the dipole operators and consider a scenario 
with small tan/3. The present bound (1) on r — /i7 translates into \D]^\ < 5 x 10"^ GeV^l In 

Fig. 13 we show contours of -D^*'"'' (left panel) and D]^^^ (right panel) in the plane {Mi,mR.^) with 
Gr = 7^/4: and rfiFi^ = 1 TcV. The D]^^'' plot refers to the choice tan (3 = 3 and n = 150 GeV. We 

^''As already emphasized, the value of the ratio (69) is quite sensitive to the parameters ^g, ^6, which depend 
on /X, the gluino mass M3 and third generation squark mass parameters. Incidentally, also affects the total 
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Figure 13: Contours of D]^"^ (left panel) and (right panel) for rriR^ = 1 TeV and Or 



IT /A. To 



determine Dl^'^^ we have also fixed tan/3 = 3 and /x = 150 GeV. 



can see that D]^"""^ and D]^^^ are comparable to or even larger than the bound in some regions 
of the parameter space, and can interfere with each other destructively or constructively. Also 
D]^'^^ may or may not exceed the bound, depending on the range of the extra parameters mf^ , 
At, Aj^^. Anyway, mutual cancellations involving all three contributions can bring the total 
dipole coefficient D]^ below the bound. This is illustrated in Fig. 14, where we plot contours 



[GeV] 




-1 
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Figure 14: Z)^ contours for tan/3 ^3,6r^ 7r/4, mij^ = 1 TeV,A^^ = and three choices of {niR^, jjL, Mi) in 
GeV: (100,150,50) (dashed), (150,200,100) (solid), (100,600,-50) (dotted). For each example the two external 
lines correspond to \DV\ = 5 x 10~® GeV~ and the middle one to Dp = 0. 



of (±5 X 10-9 GeV-^0) in the plane {A^^rhf^) for tan/3 = ^^ttir^ = 1 TeV, = and 
three choices of {fhji^, n, Mi). As a general feature the allowed regions are wider than in the 
(LFV) L examples (see Fig. 6) because the individual dipole contributions are now smaller. In 

particular, it is easier to fulfil the bound even with vanishing trilinear terms. 

In Fig. 15 we show the branching ratios of the different LFV processes in two examples, 
where fi is cither small (upper panel) or large (lower panel). In each case the dipole D]^ 



A width (through the main channel A — * bb) and hence BR{A tt) and BR{A ijlt) . 
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Figure 15: Branching ratios of LFV decays for tan/? = 3, rhu^ = 1 TeV, niR^ = 100 GeV, 9r = 7r/4, 
A^^ = and mg = 100 GeV, rhq = 200 GeV. In the upper panel, the remaining parameters are: /i = 150 GeV, 
Ml = 50 GeV, rhf^ = fhp,^ = 200 GeV. In the lower panel: n = 800 GeV, Mi = -50 GeV, = Q and 
fnjiL = "^fi,- The solid horizontal line indicates the present bound on BR(t n'^). In the upper example, 
BR{t ^ /ztt) ~ 2 X 10-", BR{t ^ /XT?) ~ 2 X IQ-^^^ BR{t jir]') ~ 4 x 10"" and BR{Z ijlt) ~ 10-^°. 
In the lower example, BR{t ^ /xtt) ~ 3 x IQ-^^^ ^^(^ ^ /xt/) ~ 4 x IQ-^^^ bR{t /xr?') ~ 10"" and 
BR{Z ^ /xt) ~ 2 X 10-". 
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Figure 16: BR{Z /it) for tan/J = 3, niR^ = 1 TeV, 0r = 7r/4 and rhfr. = 300 GeV. In the left panel: 
ffiR^ = 100 GeV and \Mi\ as shown. In the right panel: = 150 GeV. 



contribution is varied (through either Ar or m^^) and the monopole ones are essentially fixed. 
The behaviour of BR{t — > BR{t — > /.iee), BR{t — > fip) reflects the interplay between the 

dipole and monopole contributions, similarly to the (LFV)^ case. Now let us examine in more 
detail the parameter dependence of monopole contributions, considering one process at a time. 

In Fig. 16 BR{Z fir) is plotted as a function of (left panel) and in the {\Mi\,7hRg) 
plane for small fi (right panel). At variance with the (LFV)l case, the interference between 
the and contributions is now constructive. For small fi we have \Af^\ > \C^\. For large 

/i, \Af^\ is suppressed. [The A^^'^^ component grows with /i, but it remains small since tan/? is 
small. A similar comment applies to Notice that BR(Z — > /xr) docs not reach ~ 10~^, 

even for small mass parameters. This is in contrast with the (LFV) l case. 

In Fig. 17 the monopole contribution to BR{t fififi) is shown as a function of (left 
panel) and in the {\Mi\,rhji.^) plane (right panel). For \Mi\ > 50 GeV we have BR{t 
lJ,fifi)D-f=o ^ 10~^°, while values close to 10^^ can only be achieved for quite a small \Mi\. We 
recall the pure dipole contribution to BR{t — > iiiiii) cannot exceed 10^^ [eq. (38)]. 

In Fig. 18 the monopole contribution to BR(t — > //ee) is plotted in the (|/i|, mg) plane (left 
panel) and in the {\Mi\,7hji.J plane (right panel). Values O{10~^) can only be reached for 
small \Mi\. The overall BR{t fiee) can be few times 10~^ when the maximal allowed dipole 
contribution is accounted (sec also Fig. 15). 

In Fig. 19 we show the monopole contribution to BR{t — > iip) (solid) and BR{r — > /^tt) 
(dashed), either in the {\ii\,mq) plane (left panel) or in the {\Mi\,ifiR^) plane (right panel). 
In the left panel cancellation effects are visible in BR{t — > /.tTr), and are due to a destructive 
interference between A^ and box contributions. Even outside such an interference region, 
however, the r fin decay is strongly disfavoured, as its BR can hardly reach 10^^°. Regarding 
T — > lip, we see that the pure monopole contribution can exceed 10^^° for \Mi\ < 100 GeV. 
Suppose the latter contribution is 0.3 x 10""^ and the pure dipole one is 0.8 x 10~^ [maximal 
value, see eq. (39)]. Then wc have BR{t ^ /xp) ~ 2 X 10-9 or 0.4 x 10"^ depending on the 
interference sign (see also Fig. 15). As regards r fir], the BR exhibits a pattern similar to 
that of r — > fiTT, but its maximal values are even smaller. In the case of r — > fir]', the interference 
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Figure 17: IVIonopole contribution to BR{t IJ-HIJ) for tan/3 = 3, rriR^ = 1 TeV, 9r = 7r/4 and rhfi^ = 
fhjXL = 300 GeV. In the left panel: rfiR^ = 100 GeV and |Mi| as shown. In the right panel: = 150 GeV. 
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Figure 18: Monopole contribution to BR{t iJ.ee) for tan/3 = 3, niR^ — 1 TeV, 9r = 7r/4 and rhf^ = 
300 GeV. In the left panel: niR^ = 100 GeV, \Mi\ = 50 GeV. In the right panel: = 150 GeV, mg = 100 GeV. 
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Figure 19: Monopole contribution to BR{t jip) (solid lines) and BR{t /xtt) (dashed hues), for tan/3 = 3, 
ifiR^ = 1 TeV, Or = 7r/4 and rrif^ = 300 GeV. In the left panel: iriR^ = 100 GeV, |Mi| = 50 GeV. In the right 
panel: = 150 GeV,mq = 500 GeV. 
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between and box contributions tends to be constructive, and the BR could reach 10 
5.2 (LFV)ij with large tan/5 and large masses 

To discuss large (LFV)r in the large tan /3 limit we follow an approach similar to that adopted 
in Section 4.2. In this limit the Higgs-/x-T effective operators can play an important role. The 
enhanced contributions to the D'jj^''^^ dipole coefficients can be kept under control either through 
mutual cancellations (for small mass parameters) or by taking large masses. The former option 
requires some fine tuning, although this is milder than in the (LFV)l case. Among monopole 
coefficients, can be slightly larger for large tan (3 [A^*^""* is proportional to cos 2/3, which gets 
closer to -1, and A^^''^^ can play some role]. Let us consider in more detail the second possibility, 
i.e. the large mass limit (with small mA)- In this case all monopole contributions are suppressed 




Figure 20: Branching ratios of r — > /i7 (dashed lines), r — > nfifi (dotted lines) and A ^ /it (solid lines) for 
tan/3 = 40, Or = w/A, ifiR^ = 3 TeV, niR^ = 0.7 TeV and m^^ = 0.7 TeV. We have also fixed tua = 100 GeV 
and BR{A tt) = 0.1. 



and so is Z — > jjLT. The decays r nee and t ^ jip are generically dipole dominated, 
so the relations (34), (36) [and the bounds (37), (39)] hold. The dominant contributions to 
BR{t — > yU/i/i) [eq. (23)] arc induced by the dipole term Dj^ and the Higgs-mediated term 
SF]^'^, proportional to Ar [eq. (59)]. Therefore, r — >• ii/iii is correlated to both r — >• and 
A ^ fiT [11]. The processes t ^ nr],T ^ ixrj' and r — > //tt are essentially induced by the Higgs- 
exchange terms 5F^'^ 5F^'°, 5Fg''^ and so they are correlated to one another 

[eqs. (70), (71)], to r ^ 3// [eq. (69)], and to A ^ fxr. 
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In Fig. 20 we show the BRs of r — > /ij, r S/i, A — > //r in the plane (|//|, |Mi|) for 
tan/9 = 40. In the left side of the plot all decays are suppressed, in part because of internal 
cancellations both in D]^ (visible in the figure) and in Aji (not shown here, see Fig. 3 in [11]). All 
BRs increase for increasing /j, because D]^'^^ and A^^ dominate. The pattern of the contours for 
T — > /X7 is different from that of r — > 3// and A — > //r, whose rates are correlated to each other 
in the regions where r — > 3// is dominated by the Ar contribution [see eq. (60)]. In particular, 
there is a region where BR{t — > /i7) is below the present bound (1), whereas BR(A fir) 
and BR{t 3fi) can reach 10~^ and 5 x 10^^, respectively. Also notice that we have fixed 
niA = 100 GeV: BR{t S/i) is smaller for larger niA, since the Higgs-mediated contribution 
scales as l/m\. The largest values of BR{t — > 3/x) are correlated to values of BR{t fxr]) 
close to its present bound (6), and to values of BR{t — > iin) and BR(t — > iirj') of order 10~^. 

6 Conclusions 

The observation of LFV processes would be a clear evidence of physics beyond the Standard 
Model. In this paper we have focused our attention on LFV in the fi-r sector by considering 
several LFV decays, namely r — > iiX {X = 7, e~^e~ , iJ,~^ iJ,~ , p, n, rj, r]'), Z — >• /it and Higgs boson 
decays into fj,T. First we have presented a model independent treatment of such decays, then 
we have analysed them in a general unconstrained MSSM framework, without assuming specific 
relations among mass parameters. In particular, we have allowed slepton mass matrices to have 
large fi-r entries (whose origin is unspecified). This possibility is not excluded by the present 
experimental bounds [essentially that on BR{t — > /i7)], even for a relatively light superparticle 
spectrum. It is well known, instead, that the analogous situation for e-/i flavour violation is 
strongly disfavoured by the stringent bound BR{ijl — > 67) < 1.2 x 10~^^ [35]. We have imphcitly 
assumed that fi-e flavour transitions are adequately suppressed by sufficiently small e-fi and 
e-r entries in the slepton mass matrices. Issues such as the generation and the amount of LFV, 
as well as other mo del- dependent aspects, go beyond the scope of our paper and can only be 
addressed in specific frameworks (see e.g. [36, 18, 19, 37]). 

We have analysed the behaviour of the various LFV processes in the MSSM parameter space 
in terms of the relevant operator coefficients. In particular, we have disentangled the dipole 
and non-dipole contributions to the LFV decays in order to better appreciate the correlations 
among them. We have provided general formulae for the operator coefficients in Appendix 
and presented a numerical analysis with several examples in Sections 4 and 5. For the sake 
of brevity, we can try to summarize our results by schematically distinguishing a few different 
scenarios. For definiteness, consider the case of large (LFV) l- 

1. Each of the contributions to the photon dipole operator fulfils the bound dictated by 
T /i7, so no mutual cancellations among them need to be invoked. Also suppose that 
MSSM mass parameters have comparable magnitudes, with the possible exception of m^. 
For convenience, we distinguish the cases of small and large tan (3. i) For small tan (3, mass 
parameters should be larger than a few hundreds GeV to accomplish the r ^ pj bound. 
This implies that monopole contributions to all branching ratios are generically smaller 
than 10~^°. Also the Higgs-mediated contributions are small, because tan/3 is small. 
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Hence r — > /iP {P — 7r,r),r)') and Z ^ fir have BRs smaller than 10~^°. The processes 
r nee, r 3fi, t ^ fip tend to be dipole-dominated [eqs. (34), (35), (36)], and their 
BRs can be O{10^^) if BR{t fi'-f) is close to its present bound [see eqs. (37), (38), (39)]. 
ii) For large tan/9, the mass parameters have to be pushed towards the TeV range to 
accomplish the r — > //7 bound. The monopole contributions are even more suppressed 
than in the previous case. On the other hand, the Higgs-mediated contributions are 
strongly enhanced and may become comparable to or even more important than the 
dipole ones (for small m^). Hence the processes r 3fi, t ^ nrj become excellent 
signatures, together with LFV Higgs boson decays, to which they are correlated. Indeed, 
one can have BR{A — > fxr) ~ 10~^ and, for small m^, BR(t — > 3/^) close to 10~^, with 
BR{t /IT)) even larger than the latter one [and BR{t — > /xtt), BR{t /ir)') possibly 
reaching 10~^]. Regarding the other processes (and also r — > 3/x, r — > //P for large m^), 
the conclusions drawn in case i) hold again. 

2. The dipole operator satisfies the t ^ fi'j bound due to cancellations among individual 
contributions and/or because of hierarchies among mass parameters {e.g. with large /i), 
like in the examples displayed in Figs. 6 and 7. Such cases are interesting because rfiL,^ and 
M2 are not forced to be large, so monopole operators could contribute to each branching 
ratio at the level of 10~^ (or more), unless cancellations occur among different monopole 
components (which indeed could happen, as we have shown). Such contributions can 
bring BR{t fiP) (P = 71,77,77') and BR{Z fir) to values O{10~^), and an inter- 
esting interplay between dipole and monopole contributions can take place in the other 
processes. This breaks the direct correlation [eqs. (34), (35), (36)] between r — >• /X7 and 
T — > jiee, T — > 3//, T — > lip. For each of the latter processes, the BR can be smaller or 
larger than its maximal /^'''-dominated value, and can become a few xlO~^. Such values 
can be reached even in the case of strongly suppressed dipole contributions. In this ex- 
treme case r pee, r — > 3p,, t ^ p,p could be as sensitive as r ^ /i7 to super symmetric 
LFV effects. All these comments apply to the small tan j3 case. If we insist on having 
small 771^3 and M2 with large tan/3, a strong fine-tuning is required to keep the total D'^ 
below the bound (through cancellations). In this special situation, the Higgs-mediated 
contributions could push BR{t — * 3p,) and BR{t p,r]) to values O{10~^) for small uia, 
and Z-dipole contributions could help BR{Z —>■ p,T) to reach 10~^. 

The case of large (LFV)r can be summarized along similar lines, although some differences 
with respect to (LFV)^ are present (see Section 5 for details). For instance, it is easier to 
fulfil the r p,j bound on the dipole operator with small mass parameters, especially at 
small tan/3. Monopole contributions are typically smaller than in the (LFV)^ case, unless 
Ml is stretched to very small values. Barring this possibility, and assuming that r — > ^7 is 
close to its present bound, the only processes whose BRs can be 0(10-*^) at small tan/3 are 
r — * yuee, r — > 3/i, r p,p, mainly through dipole contributions. At large tan/3, also Higgs- 
mediated contributions can play an important role. If /i is larger than other mass parameters 
and niA is small, BR{t — > 3p) and BR{t — > p^rj) can reach values O(10~^). 

As a final remark, we recall that future collider experiments will give crucial information (or 
constraints) on the MSSM parameters, including those related to LFV. Indeed, supersymmetric 
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IJ,-T LFV effects can also sliow up, already at the tree level, in the decays of superparticles di- 
rectly produced at LHC or leptonic colliders [38] . This important search will be complementary 
to that of the loop-induced LFV decays discussed in our work. 
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A Operator coefficients in the MSSM 
A.l MSSM Lagrangian and conventions 

We recall here the lagrangian terms relevant for our computations, and in this way also establish 
our conventions and notations. The gauge bosons we need are the neutral ones (photon and 
Z). The covariant kinetic term for a generic Weyl fermion ijj is itjja'^D^ilj, where ip can be either 
a matter fermion (/, Z"^) or a gaugino {B'^,W'^,W^) or a higgsino {Hi, Hi , H2): and the 
covariant derivative reads as: 

^mV' = + ieQfA^ + igz{Tl - Q^sl,)Z^) ^/j, (74) 

with gz = g/cw {cw — cos OwiSw — sin^^K; where 6w is the weak mixing angle). Electric 
charge (Q^) is related to weak isospin (T^) and hypercharge {Y^) through — T^-\-Y^. Simi- 
larly, for a generic complex scalar (f) {f — /l, — Hi, i?^, Hi, H2), we have (D^0)*(£)^0). 
The gaugino and higgsino mass terms are: 

~Mi&& - ^M2#°#° - M2W^W' - fi{H^H^ - HiH+) + h.c. (75) 

Notice that our sign convention for the /i parameter is opposite to the one commonly used. The 
gaugino masses and the parameter will be taken real in the diagrammatic computations. The 
relevant fermion-sfermion-gaugino interaction terms are: 

-g'V2 (Yfff + Yf.rr - \hThI + \Hl*Hi) 

-g^ (rfff + \hI*hI - \hI*h'^ 

-g (/:/, + Hi* Hi) - g {fj^ + H^^Ht) W- + h.c. (76) 

where fu,d{fu,d) are up- or down- type components of fermion (sfermion) doublets. The /i^- 
Yukawa interactions we use are: 

-hrH^T^r - hr^H^T - KT^iHlf - Hiv^) + h.c. (77) 
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Analogous expressions hold for other Yukawa couplings (in particular, h^) . 

The (soft-breaking) scalar mass terms have the generic form — mjcl/'^P (equivalently 
denoted as —fn'j ~" ""^j I/rP)- We assume that such terms are flavour conserving, except 
for those in the smuon-stau sector, which require a more detailed description. In the left sector 
the relevant terms are: 

-Llrhl^,L, - Lirhl/L, - (Ltmi^^L. + h.c), (78) 
where = {u^, (11)^, Lr = {ur, t'l)'^ ■ In the right sector the relevant terms are: 

Finally, in the left-right sector the terms with neutral Higgs bosons are: 

where A^, A^, are mass parameters, which will be taken real in the diagrammatic compu- 
tations. 

The mass parameters in eqs. (78-79) are the entries of the following mass matrices: 

Ml = ( "^^^ "^t- ), Ml^( "^'^^^ "^^^^ ) , (81) 

where we have now assumed fh\^^^. fn']^^^ to be real. The flavour states L^, are related to 

the M.\ eigenstates L2 — {p2,(-L2V — (^Si^Ls)'^ by the relations = clL2 — slL^, = 
slL2 + clL^, where cl = cos6l,sl = sin^j;,. Analogous relations hold for the right-handed 
sleptons: = crIr^ - srIr^, tr = srIr^ + cr£r^, where £r^, £r^ are the eigenstates of Mr 
and cr — cos9r, sr = sihOr. The mixing parameters satisfy the following relations: 

~ 2 ~ 2 
SlCl = — ^^-TY- > SrCr = — ^^-ri- , (82) 

^L2 - ^L3 ^R2 - ^R3 

where and m|j^ (a = 2, 3) are the eigenvalues of M\ and M.'r, respectively. 

The mass and interaction terms presented above are sufficient to perform diagrammatic 
calculations of the operator coefficients. As anticipated in Section 3.1, electroweak breaking 
effects (Higgs insertions) will be treated at lowest order. In the next sections, for each coefficient 
we show the relevant diagrams and display the analytical results. In diagrams with gauge boson 
insertions in a fermionic line we use a shorthand notation, explained in Fig. 21. 

A. 2 Loop integrals 

The results of our diagrammatic computations will be expressed in terms of the following 
standard loop integrals 

J I 2 2\ d/c , \Af+l ^ 
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w+ _ w+ <, w+ ^ w+ w- <. w- w+ w-^w+ 



w- w+^w+ _^ w- <> w- w+ 



^1 2 ^2 _ ^2 > ^2 , ^1 ? ^1 ^2 



{H° HI) (i?? HI) {Hi Hi HI) 



Figure 21: Shorthand notation used in some diagrams. 



TT^y (fc^ — mf) . . . (/c^ — m%) 



Jo 



ds 

{s + ml) . . . (s + mjf) 



(83) 



U) 



f-1 



PC 

^0 



ds 



(s + m?) . . . (s + m%) ' 



(85) 



It is worth recalling that the expression of a diagrammatic computation in terms of these 
loop integrals is not unique, since such functions are not independent (this fact is also reflected 
in the arbitrariness in the parametrization of diagram momenta). Several identities relating 
the above functions can be derived, either by simple manipulations of the integrands {i.e. the 
propagators) or by using obvious scaling properties (see e.g. [27]). Examples of the former kind 
are 



iNiml, . . .,m%) 
JN{ml, . . .,m%) 



ml — m% 



T / 2 2\|2r/'2 2 \ 

lN-i{m^, m^_J + mj^lN{m^, m^). 



,mlf) 



(86) 

(87) 



An example of the latter kind is (for N > 2) 

d_ 



e^-%(em^,...,em^) 



0. 



By combining identities of both kinds, new ones can be obtained. For instance, from eqs. (87) 
and (88) we get 



N-2 



ET/2 22 2\ 27-/2 2 2 2\ 

i=l 

-m%lN+i{ml, . . . ,m%_^,m%,m%). (89) 
Specializing this equation to the case N — 4 and m| = m| gives a particularly useful relation 

T/2 2 2 2 2\|T/2 2 2 2 2\ o27-/2 2 2 2 2\ /nn\ 

J5(mi, mi, m^, m^) + J^^m^.m^, m^, m^, m^) = -2m^h{m-^, m^, mg, mg, mg). (90) 
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All such identities can be used to cast a given expression into a convenient form. One can 
even reduce all loop functions to a simple one^^, although the latter operation may not be the 
best choice for numerical evaluations. 

A. 3 Contributions to Af^j^ 



Hi 



Hi 



i i 



Hi Hf 



W+ IH^l W+ 



Hi Hi 



T ^ 



T 



H?HO 
W° \h?\ 



H^H^ 



T V 



HfHf 

w" \h?\ 



T V 



1^ 



HfHf 
i t 

SO \h?\ B" 



H^H^ 
t i 

so Iffpi ^0 



Figure 22: Diagrams that contribute to J^°^ (first and second rows) and (third row). 



The coefficients receive contributions of three different types: = Af^^^ + A^j^"^ + 

^L^R- The diagrams corresponding to A^j}^ are depicted in Fig. 22 and give: 



SlCl 



-(2 + 3C2^) (/xV5(M2, M|, rhlj + 2J^{Ml M|, ^x\ rhlj) 



167r2 8 

-(2 - 3c2p)Ml [i^hiMl Ml ix\ mij - h{Ml Ml /z^, 

+As20iJiM2 MMl Ml ii\ ii\ ml^ [ 



+ 



+ 



167r2 4 



C2/3 



/xV5(Ml^ M|, mij + 2 J4(Ml^ M|, ix\ ml) 

-M^M, {fi'h{MlMlf,',f,',miJ-h{MlMlf,',mlX 

^' ^" -fj:'MMlMlf,',f,',niiJ-2J4MlMlfj:',rniJ 



167r2 8 



C2/3 



15 



For example I3, which has the explicit form I^{x,y, z) = (a;y log |+j/^;log |+^;a;log ^)/[{x—y){z — y){z — x)]. 
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R 
SrCr 



{L2 ^ L3) 



(91) 



167r2 2 



C2/3 



/.V5(M2, Mf, m|J + 2J^{Ml Ml i^l rh%] 

-Ml {^I'UMl Ml rh%) - h(Ml Ml m|j); 

(i?2 ^ i?3) 



(92) 



mi 



(SO) 



r' 



HI-- 



Figure 23: Diagrams that contribute to A^^^ (left side) and A^^^ (right side). 
The diagrams corresponding to A^^^ (Fig- 23) give: 



slCl 



+ 



167r2 4 



J 4/^' [-4('^5(m2,^ 



2-2 ~2 ~2\i T/!\^2 2 2 -2 ~2 



~ 2 ~ 2 

2 ..2 ..2 ~2 



SrCr 



(1^2 i^s) 

1,2 „2 1 

If'T-^W -'-2 2r 2/t/ji^2 2-2 ~2 ~2\,t/ji^2 2 2 -2 -2 



-cl {MMI III ml^,ml^,mlj + J^{Ml /x^ ml^,mlj) 
The diagrams corresponding to ^^'^ (Fig. 23) give: 



A 



Z{c) 



"'t'^W ^ „2 ..2 



-^ 2 2 2/2T/7i,r2-2 -2 -2 -2^ 

167r2 i*^'^ (sRJ5(M2,m^^,m^^,mi^,mi^^ 



2-2 -2 -2 -2- 



(4J5(M2^ 



-2 -2 -2 -2^ 
"^K2'"^iJ2'"^L3'"^L3> 



(93) 



(94) 



+4J5(M|,m' 

+(^RMMlm%^,mji^,ml^,mlj') - (4-4) (4■^5(^2^ "^L' "^ia' ^L) 
+c\j^{Ml m%,m^R^,ml^,mlj) 
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+ 



U2 2 



1 2 2 \ 2 ( 2 J / T{/r2 ~ 2 ~ 2 ~ 2 ~ 2 \ 



167r2 4 

+clMMl m^3, m^3, mL, mij) - 4 (s^slMil : 

i2t/7i^2~2 ~2 ~2 ~2\\ / 2 2\(2i(-\jr2~2 ~ 2 ~ 2 ~ 2 \ 



~ 2 ~ 2 ~ 2 ~ 2 ^ 



(95) 



A 



Z{c) 
R 



2„2 



Kis 



SrCr 



T^W 2 2\ 2(2 1 /T\/r2 ~2 ~2 ~2 ~2 \ 
i2t/7i^2~2 ~2 ~2 ~2\\| 2/2t/7i^2~2 ~2 ~2 ~2\ 

+Cl M^i, rriL^ , 771^3 , TTiR^ ,mR^)j+CR [s^ J5 (M^ , m^^ , m^^ , , ) 



+Cl^5(Mi , m^g, m^g, ^i?3 
A. 4 Contributions to C2 ij 

The diagrams corresponding to CI^r are depicted in Fig. 24 and give: 

,2 1 



Cl_ 
SlCl 



9 



167r2 12 

-4i^5(M2, Ml mlj + mlh{Ml Ml Ml Ml mij 



+ 



CI 

SrCr 



167r2 12 

9" Ir 



167r2 3 



-K5{Mlm%,m%,m%^,m%^)] - {R2 R3 



(96) 



(97) 

(98) 



>7 



(BO) 



5° 



n ^7 
Figure 24: Dia grams that contribute to C2 ^• 



A. 5 Contributions to 

The diagrams corresponding to are depicted in Fig. 25 and give: 
0^ 1 



167r2 24 



( 1 + 2s^ ) (M^^ mi^ , mi^ , mi^ , mi J 
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SrCr 



+ 



-4c2^ {2K,(Ml Ml Ml Ml mij - SM^^Ml Ml Ml M|, Tf^ij) 
-1 + 2slr)MMlml^,ml^,ml^,7filJ - {L^ ^ L3) 



9" 1 
167r2 24 



■ ^sl,K5{Mlm%,m%,rh%,m%) - {R2 R3) 



(99) 
(100) 



< w+ & 



Figure 25: Diag rams that contribute to C^^- 



A. 6 Contributions to D\ji 



The coefficients DI^r_ receive contributions of tfiree different types: -D2,i? = ^1[r'^^l%'^^l%- 
The diagrams corresponding to D]}'^ are depicted in Fig. 26 and give: 



^^^^ - ^ mi^J5(M2^mi^,mi^,mi^,miJ-2M2V5(M2^M2^M2^M2^miJ 



slCl 167r^ 4 



+ I|^■i^i2>^5(Ml^mi,,7f^i,,7f^i,,miJ-(L2^L3) (101) 



^^^^ - ^L^5(M2,m|^,m|^,m|^,m|J-(i?2^i?3) (102) 



>7 



C V 



iR. y 



7 



Figure 26: Diagrams that contribute to -D^^^^. 
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The diagrams corresponding to are depicted in Fig. 27 and give 



,7(6) 



167r2 2 



41 



D 



,7(6) 
R 



SrCr 



+ 



+M2//tan/3 (mlj^{Ml mi^, mi^, mij + 2J^{Ml M|, ii\ mij 
+2 J5(M2^ M|, mij + 2 J5(M|, /x^ //^ /x^ mij 



9' 



12 



m 



L2 



167r2 2 

(i^2 ^ i^3) 



(103) 



(104) 




7, 



W+ 



Hlk (SO) 



-«2t (SO) 





?1 






— * \ * * — • — * — 








Hi \ 





7< 



Hj;)H+\ W-W+ 



7 



7 



I7 



I7 



Figure 27: Diagrams that contribute to D]}''^. 
The diagrams corresponding to -D^*^^^ are depicted in Fig. 28 and give: 



D 



.7(c) 



167r2 ^ 



{ [sLCL {s\[Ar + tan/3] + srCrA"^^) + 44^^^ /5(Ml^ mi^, m|J 

- [siCi (4[A. + tan/3] + Si^c^Aj,) - 44^^ h{Ml ml^,mlj 
+ [slCl (4[^r + tan/3] - SkC^AJ,) + 44^Jr] ^5(M2, mi^,mlj 

- SLCL (4[^. + /Xtan/3] - SrCrA^,) - 44<r] /5(M2,M2,M2, 771^3, 



(105) 



D 



7(c) 



167r2 



h{Mt,Mf,Mt,mi^,mj,^ 



+ 



SRCR [sl[Ar + //tan/3] + SiC^Aj^) - 44 Aj^] /5(Ml^ M^, M^, ml^,m%) 
srCr {cl[A, + //tan/3] - s^Ci^lJ,) + 44<r] ^5(Ml^ 7h%,7h%) 
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SRCR [cl[Ar + At tan/3] - slClAJ^) - 44 Aj^] h(Ml rfil^, ml^)} 



(106) 



BO 



BO 



BO 



BO 



BO BO bo bo 



Figure 28: Diagrams that contribute to -D^ (first row) and D^^'^' (second row). 



c 



A. 7 Contributions to Df 



L,R 



The leading contributions to D^r are of two types: D^^ r — D^r + D^r. The diagrams 
corresponding to D'^^r are depicted in Fig. 29 and give: 



D 



Z{b) 



SlCl 



D 



m 

R 



SrCr 



+ 



- M2// tan/3 -2(1 + 2si^)mijr,{M^, rhi^,rhi^,rhi^) 



167r2 8 

+2(1 - Asl,)MMl //^ //^ //^ mlj + Sc^^MMI Ml Ml //^ mij 
+ (5 - 85^^) J5 (Ma^ , Ml ,fi\fi\ ml^ ] 



9 



/2 



1 



167r2 [-2(l-2s^)mi^/5(M2,/x2,mi^,mi^,miJ 

+2 J5(Ml^ /x^ /x^ mij + J5(Ml^ //^ mij] - (L2 
5'^ 1 



-Ml// tan /3 -As^rfiRj^iM^, ^ ,mR^,rhR^,rhRj 



167r2 4 

+2 j5(m2, /x^ /x^ 7^4) + j,{Mi Ml ^\ fi', 7^4; 



iR2 



L3) (107) 



Z(c) 

The diagrams corresponding to r are depicted in Fig. 30 and give 



D 



Z{c) 



9 



1 



SlCl 



167r2 4 



-Mi/xtan/3 (1 - 2s^) [s'rMMI Mlrhi^,7hi^,m'R^ 



+clMMl Ml rhl„rhl„m%)) - 2s^ {s^MMI Ml m^^, 
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;^2t 


{&>) 
















~ *- ' 















tJ-O : 



fx 



H^H+ \ w-W+ 



rrO : 






^0 B° 







T ^ 



Figure 29: Diagrams that contribute to D^j}^^- 



D 



Z{c) 
R 



SrCr 



+cj,MM^, Ml mi^, m^3, m'^^)) - (L^ ^ L3) 



,/2 



(109) 



167r2 4 

+clMMlMlml,m%,m%j)] - {R^ ^ iis) (HO) 



^z \m 



J^2 : 



B" 



^z \m H^. z^ 



Figure 30: Diagrams that contribute to D 



Z{c) 
L,R ■ 



A. 8 Contributions to B 



fL,R 

L,R 



Contributions to B)^'^ arise from the diagrams in Fig. 31. 
For f i^,t: 



1 



(1 + 4 5^3^,-1) MMi, Mi,rhl^,rh\) 



slCl 16n^ 16 

-2 (1 + 45^3 ^^]M■^h{MlMlml^,mjJ 



+ 



ivl"^ ^ / fl/f2 fl#2^2 ^2 
y y T-i3 \ T /'7i/r2 n^2 ^2 ^2 



■ i^i-TfA) [HMl Ml ml, rhjj + 2M,MMMt, Ml rhl,mjj 
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1 »- 1 

W+ W+' 






»- 1 

■W+ w+' 

w- W-- 
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— »■ — 1 
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». 1 
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1 K 1 


1^ 
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K 1 
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» 1 
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h 


f 


r 


>- 

Ir 
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Figure 31: Diagrams that contribute to S^^'^ (first and second rows) and Bj^ "^ (third row). In the case of 
B'l'^R, i-e- f = = f = the following replacements are implied: fl v^Jl ^ ^i^, /fl ^ ^jj^- 



+ 



B 



Ir 



SlCl 



SrCr 



167r2 4 

o"^ 1 
167r2 4 •'^^ 



J4(M/, M/, mi,, m)J + 2Mth{Ml, Mf, ml,mjj 



{L2 

,/4 



9 y2 



167r2 



B 



fR 

R 



/4 



SrCr 



167r2^ ^/fl^ 

{R2 ^ i?3) 



For / = A*: 



-4 (5J4(M2^M2^mi,,miJ +2M2/4(M2,M2,mi,,mi, 



(111) 



(112) 



(113) 



(114) 



slCl IBtt^ 16 

+4 {5MMlMlml^,ml^) + 2Mlh{MlMlml,mlj) 

+{cl - si) {5MMl Ml ml,ml) + 2Mlh{Ml M|, ml^.ml;)) 
+ t^-l \-cl {MMl Ml mi,, mij + 2M,M2h{Ml M|, mi„ mij) 



167r2 8 

+4 {MMlMi,ml^,7hl^) + 2M^M2h{MlMlml,ml)) 
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slCl 



SrCr 



B\ 



+(cl 



^ • hi (J4(M,^M,^mi^,<) + 2M,2/4(M^M^<,mij) 
+4 (J4(Ml^Ml^mi3,mi3) + 2M2/4(M2,M2,mi3,mi3)) 
+(4 - 4) (j4(M^M^<,mi3) + 2M^74(M^M^mi^,mi3)); 

= ^ • ^ [4 (^4(M2, M^, mi^, m^J + 2M2/4(M2, M^, mi^, m?, 1^ 

9/t/i.^9 ,.^9 ~9 ~9\ .^,.^9^/1.^9 ,.^9 ~9 ~9\\" 



=3))] 



iDTT' 4 L - V - 

+4 [UMl Ml ml,mlj + 2M^h{Ml Ml ml^,m% 

^ ^4 [j,{MlMlml,m%) + 2M',h{MlMlml, 

! 71^2 .~2 .~2 \ , r,nj-2T/nj-2 71 ^2 .~2 2 



- (i^2- 

167r2 4 L - V - - ^ - - 

+4 {MMlMlml^,mlJ + 2Mih{MlMlrhl,rhl^ 

+ 2M^h{MlMlm% 
^i,m^3,m^3)) 



-4 



(j4(M2,M2,m2^^,m^^ 

(J4(M2, Ml ml^,mlj + 2M^h{Ml Ml m\^,m 
+(4 - 4) (^4(Ml^ Ml m\^,m\^) + 2M^h{Ml Ml 



+4 



(115) 



(116) 



(117) 



rh%,mlj)] (118) 



A. 9 Contributions to Al^r 

The leading contributions to are of two types: A^^^ = A^^''fj + A^^^^^. The related diagrams 
are depicted in Fig. 32 and give: 



SlCl 



- (L^^Ls) 



4^^ 



9' 



12 



SrCr 



167r2 



piMih{Mt,ixlmi,J - (i?2 ^ i?3) 



(119) 
(120) 



A 



R 



SrCr 



J2 



^/.Mi [-4^3(M,^mj,^,miJ - c'^h{Mlm'^^,mD\ - {U - ^3) (121) 



-,'2 



167r2' 



/.Ml -sih{Mlml,mj,^)-cih{Mlml,mi^) - (R, ^ R,) (122) 



A. 10 Contributions to muon g — 2 

The MSSM contributions to the muon anomalous magnetic moment are described by diagrams 
analogous to those for D'^ (Figs. 26, 27, 28), with the replacements t(t'^) — > //(//'^). The result 



46 



2i 



SO 



can 



a 



Figure 32: Diagrams that contribute to (first row) and A]^'^ (second row). 

1 be decomposed into five terms, i.e. a^^^^ — a^^L + a!"^R + + o!^^\ + a^^^ given as follows: 



2m2 



o2 2 

167r2 4 L ^2 2, L 



■2M|J5(M2,M2, 



Ml Ml 



-2 ~ 2 



+ 4 



/2 



167r2 



(a) 

2ml 
2m2 



r2 

^L~2T(')i/r2~2 ~2 ~2 ~2\|/ 



„2 r 
Sl, -^2 



r7'2 

y 2~2t/7i^2~2 ~2 ~2 ~2\|/2 2d 



(123) 
(124) 



167r2 ' 2 
+M2Ai tan 13 



-rhl^MMl /i^ mi^, mi^, mij - 2mi^75(M|, 

mi^/5(M2^//^m2 



2 ~ 2 
// ,m 



+ 7^ 



^2', A*', ^L) + 2 J5(M|, /x^, f/, ml^ ^ ^ 

/2 2 

^ y ["^5 (^1 ' ^^2 ' <2 , mi J - Mi/i tan /3 /5 (M^ , /x^ , 
(rH .s? . Ln L3) 

J5(M,^ - 



2 ~ 2 ~ 2 \ 
'L2'"^L2'"^i:2) 

r2 2 2 ~ 2 \ 



"'fiR _ 9 2-2 

2m2 - IGtt^'^^'"^^ 

+ (4 ^ i?2 

2m2 167r2 ^ 



2 ~ 2 ~ 2 

^i?2'"^R2'"^ii2 



■i?3 



J +MlA^tan/375(Ml^ 



(125) 



2-2 - 2 - 2 > 



(126) 



(c%slClA^^ + srCrc\A^^^ + srCrSlCl[A^ + /xtan/3]) 

h{MlMlMlf~'' -2 ^ 



+c\c\{A^ + /xtan/?) 



- 2 - 2 \ 
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+ 



— {-CrSlClA^^ + SrCbsIA^^ - SRCRSLCL[Ar 

+44 (^M + A* tan P)] h {Ml , Ml , Ml , ml^ , m% ) 

— [s^slClA^^ - SrCrcIA^^ - 



- /itan 



ir - SRCRSLCiiAr + /xtan/?]) 
+44 {A^ + tan P)]h{Ml, Ml , Ml , ml^ , m% ) 

777- / 

+ — (sj^SLCLA^^ - SrCrsIA^;^ + SrCrSlClIAt- + tail 

+44(A^ + A* tan/3)] h{Ml, Ml, Ml, rhl^,mlj} 



(127) 
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